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ABSTRACT 

We examine the "puffed-up inner disk" model (Dullemond, Dominik & Natta 2001), proposed for 
explaining the near-IR excess radiation from Herbig Ae/Be stars. Detailed model computations show 
that the observed near-IR excess requires more hot dust than is contained in the puffed-up disk rim. 
The rim can produce the observed near-IR excess only if its dust has perfectly gray opacity, but 
such dust is in conflict with the observed 10/^m spectral feature. We find that a compact (~ 10 AU) 
tenuous (ry < 0.4) dusty halo around the disk inner regions contains enough dust to readily explain 
the observations. Furthermore, this model also resolves the puzzling relationship noted by Monnier 
and Millan-Gabet (2002) between luminosity and the interferometric inner radii of disks. 
Subject headings: accretion, accretion disks — circumstellar matter — dust, extinction — stars: 
pre-main-sequence 



1. INTRODUCTION 

Processes in the immediate vicinity of young pre-main- 
sequence stars influence the initial stellar evolution and 
the formation of terrestrial planets. Since small scales 
of several AU around a star are difficult to resolve, we 
still lack a clear understanding of processes such as disk 
accretion, launching of bipolar jets and winds, and dy- 
namics and reprocessing of dust in the inner hot disk re- 
gions. The dust geometry is one of the basic ingredients 
needed for constraining theoretical models of these pro- 
cesses. Traditionally, this geometry has been deduced 
from the spectral energy distribution (SED), which is 
dominated at infrared wavelengths by dust emission. 

A widely popular geometrical description is 
the two-layered flared disk model developed by 
iChiang fc Goldreichl lfl997|) (CG hereafter). The model 
gives a simple method for estimating the flux from the 
optically-thin surface layer of an optically-thick disk 
directly exposed to the stellar radiation, and from the 
colder disk interior heated by the warmer surface. The 
simplicity of the method, together with evidence for 
the existence of disks based on radio imaging, made 
this model a dominant description of T Tau and Herbig 
Ae/Be stars (intermediate-mass, 1.5< M*/M Q <10, 
counterparts of T Tau; Haebes hereafter) . 

Although the CG model successfully explains the 
observed SEDs, advances in imaging techniques re- 
vealed shortcomings of this model. Analyzing im- 
ages at scattering and d ust emission wavelengths, 
IMiroshnichenko et al.l l)1999f) concluded that disks alone 
cannot explain the imaging observations, at least not for 



some Haebes. Instead they modeled the SED with an 
optically thin halo surrounding an optically thick disk, 
and emphasized that only multi- wavelength imaging can 
distinguish between this and the CG model. Subsequent 
detailed modeling of imaging data in numerous systems 
revealed the existence of dusty halos around the putative 
flared disks llGo mez fc D'Alessiol 120001: iSchneider et alJ 
12003 iStanelfeldt et alJ l2003fl ~The inadequacy of the 
SED as the sole analysis tool in d etermining the geom- 
etry was further demonstrated by iVinkovic et all l)2003fl 
(V03 hereafter). They showed that the mathematical ex- 
pression for the SED calculation in the CG model can be 
transformed into that for the halo-embedded-disk and 
vice versa. This has far reaching consequences for all 
studies based solely on SEDs. If not supported by imag- 
ing at various wavelengths, SED models can lead to erro- 
neous conclusions about the spatial distribution of dust. 

The disk inner region in Haebes (within ^10 AU from 
the star) proved to be more complicated than the orig- 
inal CG model. Thermal dust emission from this re- 
gion peaks at short wavelengths, creating a near IR 
bump (l/j,m < A < 8/xm ) in the SED o f many Haebes 
(Hillenbra nd et alJll99l . IChiang et all l|200ll) noticed 
that the CG model did not produce enough near IR flux 
to explain the bump. This implies that the disk flaring, 
which increases the emitting volume of the optically thin 
disk surface, is too small at the inner radii. Since the disk 
geometry is constrained by vertical hydrostatic equilib- 
rium, an additional hot dust component is required for 
explaining the near IR bump. To s olve this problem, 
iDullemond. Dominik. fc Nattal 1)200 lft (DDN hereafter) 
proposed to modify the CG geometry without introduc- 
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ing an additional component. They noted that the disk 
vertical height is increased (puffed up) at its inner rim 
because there the disk interior is directly exposed to the 
stellar radiation and hotter than in the CG model at the 
same radius. The rim is the hottest region of the disk 
and with its increased size it is possible to boost the near 
IR flux. This puffing of the rim is e quivalent to the disk 
extra flaring that was identified by iChiang et al.1 (|2001f) 
as missing in the CG model. 

Evidence in support of the DDN model was gar- 
nered from SEP mod eling of a large sample of Haebes 
( Dominik et al]l2003ft and recently also of T Tau stars 
i Muzerolle et al.l 1200.*$) . Still, the inner disk geome- 
try remains controversial. Recent advances in near 
IR interferometry provide imaging data of this region, 
and the first results from a large sample of Haebes 
show that man y of these objects appear close to circu- 
lar symmetr y (Millan- Gabet. Schloerb. fc Traubl I200H 
lEisner et al.ll2004j) . This" is an unusual result if disk in- 
clinations are random. It also creates a new set of prob- 
lems when interpreted as almost face-on disks because 
that often conflicts with outer disk inclinations derived 
from other imaging observations (HST, radio). This is 
difficult to accommodate in disk-only models, but is eas- 
ily explained by halo-embedded disks (V03). 

In this paper we reexamine the DDN model and the 
theoretical approach behind it, and identify some unre- 
solved issues in its description of the rim emission. We 
employ exact radiative transfer calculations of the rim's 
brightness and show that the concept of puffed up rim 
requires some fine tuning of the model parameters in or- 
der to produce enough flux to explain the observations 
(e.g. the dust must be perfectly gray). 

Various independent observations indicate the exis- 
tence of compact halos (~10 AU) around the disk inner 
regions (see V03 and references within) , and we find that 
such halos readily explain also the observed near-IR ex- 
cess. Furthermore, the halos also resolve the puzzlin g re- 
lationship noted bv lMonnier fc Millan-Gabetl J2002T) be- 
tween luminosity and the interferometric inner radii of 
disks. 

2. EMISSION FROM THE INNER WALL 

A distinct feature of the near IR bump is its anoma- 
lously high flux F\ as compared with the stellar emis- 
sion. To quantify this effect we introduce the flux ratio 
^/jm/fi/im as a measure of the strength of the near IR 
bump; this ratio increases when the prominence of the 
bump becomes larger. These wavelengths are chosen be- 
cause the 2fim flux is dominated by the dust, while the 
1/xm flux is dominated by the star. 

Figure Q summarizes the observed values of 
the -F2/im/-Fi/im flux ratio for a sample of well- 
observed stars, with data c ompile d from the fol- 
lowing references: lLow I (1197(1: [Gillett fc Stem] 
( 1971 : JStrqm et all (U972I) : lAllenl (I1973D: iCohenl 
(Il973albllcl): IGlass fc Periston! 111971: IGoherJ (119751): 
ICohen fc Schwartz! lll97fift:lKolotilov et alJ (ll977l):lCohen 
(1980); Bouchet fc Swings] (1££2J); Iborenzetti et~ 
( 1 1983ft : iTiin A Diie et alJ (119841): Kilkenny et al. 
(1198!: iThe et a.lJ fl98l: IQlnon et al] (l198fl): 
IBerrilli et alJ (I1987I): [Strom et all (I1989D: IHu et all 
(11989ft: lLawrence et all (I1990D: iFououe et alJ dl992fl: 
IBerrilli et alJ l)1992ft : iHutchinson et all (|1994j) : iLi et all 



(11994ft: iPrusti et all (I1994H : iSvlvester et~all \199& : 
Garcia-Lario eHU IMalfait et alJ (TS: 

Herbst fc Shevchenkd (11999ft: Ivan den Ancker et alJ 
( 2000): Ide Winter et all (|2001(l and A. S. Mirosh- 
nichenko (2005, private communication). 

Dust extinction at lfim is larger than at 2/im and could 
enhance the observed strength of the near IR bump by 
^20% for Ay — 1, therefore only objects with Ay < 1 
were considered. Since the reddening correction is negli- 
gible, the uncorrected data displayed in the figure repre- 
sent the true range of near IR bump strength in Herbig 
Ae stars. The underlying stars of all objects have tem- 
peratures of about 10,000 K, which gives -F2/um/ Fi/im 
= 0.09. Yet in all objects this ratio exceeds 0.25, re- 
flecting a large NIR excess from hot dust emission (Hil- 
lenbrand et al. 1992). The luminosity of each object is 
displayed together with its name in figure^and it ranges 
from ~5L Q to ~80I/©. Luminosity does not show any 
correlation with the near IR bump strength, reaffirming 
our conclusion that these data can be used as a general 
description of the near IR bump strength in Herbig Ae 
stars. 

2.1. General Description of the Rim Emission 

At the inner rim, gas that is typically part of the disk 
cold interior becomes directly exposed to the stellar ra- 
diation and expands to higher scale heights. According 
to DDN, emission from such a puffed-up rim can explain 
the near IR bump in the spectrum of Herbig Ae/Be stars. 
The rim geometry is sketched in figure[21 The rim is mod- 
eled as a cylinder of radius R r im and height 2H r i rn cen- 
tered on the star. The basic assumption of this model is 
that the rim is optically thick in the near IR and shorter 
wavelengths. This maximizes the rim energy output. 

The original DDN model 

(jDullemond. Dominik. fc Nattal 12001ft successfully 
explained the data, but it was based on an approximate 
treatment of the rim height and emission. More realistic 
models were calculated by the authors o f the DDN 
model in their subsequent work. iDullemondl ({2002) used 
a 2D radiative transfer model for gray dust combined 
with the hydrostatic equilibrium. The obtained near 
IR bump strength is shown in figure ^ (solid line) . The 
maximum strength is still too low to explain all the 
data, but it can accommodate the majority of observed 
near IR bump strengths. 

A dramatic reduction of the DDN model efficacy hap- 
pens w hen a mixture of small and b ig grains is intro- 
duced. iDullemond fc Dominikl (|2003)(DD04 hereafter) 
combined 2mm (big, gray grains) and 0.1/^m (small) 
grains in various ratios and performed 2D radiative trans- 
fer calculations coupled to the equation of vertical hydro- 
statics and dust settling. The model fails to explain the 
data even when 99.999% of the dust mass is in big grains 
(see figure ^1 . The behavior of this result is unexpected; 
decrease in the small grain population leaves more gray 
dust grains in the mix, which should move the whole so - 
lution closer to the gray dust result of lDullemondl {2002) . 

A closer inspection of obtained results shows that the 
temperature of big grains in such a multi-grain mixture 
is lower than in the pure gray model. In a mixture, both 
small and big grains absorb a fraction of the local en- 
ergy density and participate in providing the local dif- 
fuse heating. But, as shown in the next section (see also 
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Fig. 1. — Strength of near IR bump in Herbig Ac stars. The "naked star" arrow on the right axis marks the strength for a 10,000 K black- 
body spectrum. Other arrows indicate the data (see J2]for refer ences) for sy stems with unknown disk inclination angles. Thes e angles were 
estimated for three star s as in dic ated: HD 163296 — 1: IGradv et aT]|200(J). 2: [ Mannings fc Sargenti il99l): HP 100546 — 1: IGradv et alj 
i2(fuT):IAugereau et ap EOOll). 2: ILiu et all <2003T); AB Aur — 1: [Eisner et alJ 12004) andlEisner et all 12003) . 2: ISemenov et all 12004). 3: 
IFukagawa et all 120oJT 4: IGradv et all 119991) . 5: ILiu et"aT] 12003) . 6: IMannings fc Sargent! 119971) . 7: ICorder et all 120051) . 8: IPietu et all 
12005). Stellar luminosities in units oi Lp, are in dicated together with the stellar name. The solid line is a 2P radiative transfer model 
of puffed up gray dust wall by|P ullcmond (2002). The filled circles with letters (realistic models, inclination i = 45° ) are 2P radiative 
transfer models by Pullemond & Pominik (2004) that include a mixture of big and small dust grains. The rim height in all these models 
is calculated directly from the equation of vertical hydrostatics. While the gray dust model can explain the majority of data, the realistic 
models of puffed up rim fail to explain the data even when 99.999% of dust mass is in big (gray) grains. 




Fig . 2. — Sketch of the puffed up inner disk wall (see M2.1l and 
i2"2lfor details). 



equation IA8I in the appendix) , small grains are a very 
inefficient source of local diffuse heating, resulting in less 



efficient heating of the big grains then in the pure gray 
model. With such a temperature decrease, the vertical 
hydrostatic equilibrium cannot produce disk puffing com- 
parable to the gray model. While in the gray model the 
puffed-up disk rim height is close to H r i m / R r i m = 0.2, 
multi-grain models have only H rim / R rim < 0.15. Since 
the observed rim emission scales with rim height, this is 
the major reason behind the failure of the multi-grain 
models to explain the data. The presented model with 
the lowest fraction of small grains yields the largest dis- 
crepancy because it suffers the largest reduction in the 
small grain contribution to gas heating and rim puffing 
while still having enough small grains to suppress heat- 
ing of the big grains. In quantifying this effect, the mass 
ratio between big and small grains that DD04 used is not 
the most illustrative choice. A more appropriate quan- 
tity would be the "equivalent" grain size of the grain 
mixture. In the case of two grain populations with sizes 
abig and a sma ii and fractional number densities X^ig and 
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Xsmaii (such that X bi g + X s 
size obeys 



1), the average grain 



(a 2 ) — Xb ig al ig + X sma ua 2 srnaU . 



(1) 



The number fractions can be deduced from the reported 
parameters of the DD04 models: fixed inner and outer 
disk radius, fixed total disk mass, and total dust mass 
in big and small grains. The model with 99.999% of the 
dust mass in big grains has only Xug — 1-25 x 10 -8 , and 



0.25/im. This grain size 



X smaU ~ 1, yielding (a 2 ) 1 / 2 
is too small to be considered equivalent to the gray dust 
model. 

It is important to note that this equivalent grain is just 
an indicator of the overall solution and cannot be used as 
a general replacement (average or synthetic) grain f or the 
radiat ive transfer calculation. As already shown bv lWoli 
( 2003) , the approximation of an averaged single grain 
as a replacement for a dust mixture breaks down at the 
surface of a dust cloud (or in this case the rim surface). 
A more detailed study of multi-grain disk models will 
be presented in a separate publication, while in the next 
sections we explore the limits of possible DDN model 
applicability in the context of single dust grains. 

2.2. Approximate Solution for the Rim Emission 

Denote by i?* and T* the stellar radius and tempera- 
ture, respectively. At distance d and direction i where 
the star is free of rim obscuration (see figureEJ), the over- 
all observed flux at wavelength A is (R*/d) 2 irB\(T r ) + 
F\ lm {})- If I" m (i) is the rim surface brightness in the 
observer's direction then 



A 



(*) = "*(*) ^H rim R rm sinz. (2) 



Here the cylindrical visible surface is replaced with a flat 
rectangle. This approximation maximizes the flux since 
curvature decreases the projected area of portions of the 
visible surface, reducing the observed flux. We also as- 
sume that the stellar illumination is perpendicular to all 
portions of the rim. This too maximizes the observed 
flux. 

The observed rim flux in equation [5] is determined by 
the rim height, surface brightness and radius. Our 2D 
radiative transfer calculations described in W2.'6\ confirm 
that the rim emission is indeed proportional to the rim 
height, therefore we maximize the rim emission in this 
study by using H r i m = 0.2R r i m , the maximum height 
allowed before the rim starts to shadow large portions of 
the disk (DDN). The solution for any other rim height 
can be derived from our models by a simple scaling of 
the rim emission. 

The surface brightness of a gray dust rim can be ap- 
proximated with B\(T r i m ), where T r i m is the dust subli- 
mation temperature. The description of a non-gray sur- 
face must take into account the spectral variation of op- 
tical depth of th e emitting optically- t hin s urface layer. 
This was done bv lChiang fc Goldrcichl l)1997ft . According 
to their model, the surface layer vertical optical thick- 
ness is unity at visual (a characteristic wavelength of 
the stellar radiation absorption) Ty = 1, therefore at 
all other wavelengths it is t\ = a^ bs /a v bs = q\. The rim 
emits at near IR where q\ < 1 (the dust NIR opacity is 
smaller than at visual), thus the surface layer is optically 



thin at these wavelengths and its emission is reduced ac- 
cordingly. Therefore, the rim surface brightness becomes 
~ q\ B\(T r i m )/ sin i and the observed rim flux is 



Fx 



f/ 2 



B\ (T r , 



II,. 



R, 



q\ non-gray dust 
sin i gray dust 



( 3 ) 

This result shows that a non-gray rim creates a smaller 
IR excess than a gray opacity rim. In addition, non- 
gray opacity removes the angle dependence from the rim 
emission (we expect this approximation to break down 
at very small inclination angles where q\ ~ sini). 

The rim radiu s is derived from radia tive equilibrium, 
which gives (e.g. Ilvezic fc Elitzurlli997t) 



T, 

Trir. 



g(r«) 



1M1 + 



H r 



"I 1/2 



(4) 

Here <r(T) is the Planck average of a^ bs at temperature T, 
ip describes the correction for diffuse heating from the rim 
interior 1 and 1 + H r i m / R r i m is a correction (described 
by DDN) for self-irradiation from the other side of the 
rim. In appendix E] we derive an approximate solution 
which shows that gray dust, with a(T*)/a(T r i m ) = 1, has 
ip ~ 4 and that non-gray dust, with cr(T*) / a{T r i m ) > 1, 
has ip ~ 1. Note that for gray dust this makes equation^] 
identical to the original DDN expression (their equation 
14). The approximate near-IR bump strength is given in 
equation IA1 01 yielding 



F", 



2/xm 



' 0.23 



non-gray dust 



i%tm [ 0.09 + 0.52 sin i gray dust 



(5) 



for T*=10,000 K, T„ m =l,500 K and H„ m / R„ m =0.2. 
Comparison of this result with the data in figure^shows 
that the NIR bump of non-gray dust is too small to ex- 
plain the observations. Therefore, interpretation of the 
NIR bump in Hcrbig Ae stars with inner disk puffing 
places a strong constraint on dust evolution in this re- 
gion. The dust must grow to a size greatly exceeding 
the initial interstellar size distribution, and small grains 
must be depleted to such a large extent that the inner 
disk opacity can be considered gray. In the next subsec- 
tion we employ exact 2D radiative transfer code to obtain 
accurate values for ip and place more precise constrains 
on the DDN model. 

2.3. Exact Models for Single-Size Grains 

To examine the validity of conclusions based on our ap- 
proximate solution we performed full 2D radiative trans- 
fer calculations for an optically thick torus centered on 
a 10,000 K star. The torus cross section is sketched in 
figure |2 it is a square with side-length of 2H rim , where 



//,. 



0.2i? r 



This configuration is the same as de- 



scribed by DDN, where the puffed-up disk rim is a cylin- 
drical surface directly exposed to stellar radiation, while 
the rest of the inner disk is in its shadow. The dust 
has sublimation temperature T su b =1,500 K and con- 
stant density everywhere in the torus, with horizontal 
and vertical optical depths Ty = 10, 000 in visual. Dif- 
ferent density structures do not change our results as 



1 Note that Ilvezic k. Elitzuil 119971) used * : 
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Fig. 3. — Same as figure but for single-size grain models with H r j m = 0.2R r i m . The shaded areas are regions without flux 
contribution from the rim, because of either the absence of dust or rim self-shadowing. The dashed line is the o rigin al DDN model 
(Dullcmond, Dominik, & Natta 2001); all other lines show the results of exact 2D radiative transfer calculations (see H2.3I . 



long as the ry = 1 layer on the illuminated surface is 
geometrically much smaller than H r i m . 

Radiative transfer modeling was conducted with our 
code LELUYA (http:/ / www deluya.org[ | that works with 
axially symmetric dust configurations, ft solves the inte- 
gral equation of the formal solution of radiative transfer 
including dust scattering, absorption and thermal emis- 
sion. The solution is based on a long-characteristics ap- 
proach to the direct met hod of solving t he matrix version 
of the integral equation ijKurucz |ll969fl . 

The results are shown in figure [3] together with the 
original DDN solution (dashed line). Our 2D model re- 
sults for gra y dust witho u t scat tering (solid line) are very 
close to the iDullemondl l)2002f) results, shown in figure 
^ which also included vertical hydrostatics equilibrium. 
This model has ip =4. Its rim radius (49i?*) and flux 
are essentially the same as the original DDN model, con- 
firming that a puffed-up rim of gray dust is capable of 
explaining the near 1R bump. 

Since realistic dust is not gray at all wavelengths, we 
calculated models for silicate dust with diffe rent grain 
radii, employing optical constants from iDorschner et al.l 
(1995) (x — 0.4 olivine). Figure |3| shows results for three 



representative grain radii, with the corresponding rim 
properties summarized in table Our m odel results 
for 0.1 /im grains are almost identical to the Dul lemondl 
(2002) results for purely small grains (see figure QJ. As 
is evident from figure the model can explain the data 
when the grain radii are 2/im and 0.5/zm, but it starts to 
fail as a general explanation of the near IR bump when 
the grain radius drops below ~ 0.1/im. 

A decrease in grain size has two opposing effects on 
the rim flux. On one hand, the ratio a{T*) / /j{T su i,) is 
increasing, leading to a larger rim radius and emitting 
area and thus enhancing the rim emission. On the other, 
the rim surface brightness is declining because q\ is de- 
creasing, reducing the rim emission. The net result is 
that maximum rim emission occurs at grain radius of 
~0.5/im, which, as is evident from figure 01 corresponds 
to the transition between gray and non-gray opacity in 
the near IR region. This is predominantly a grain size 
effect; the dust chemistry introduces only second order 
corrections. 

When the grain radius drops below 0.5/Ltm the dust 
opacity becomes non-gray in the near IR and the puffed- 
up rim model begins to fail. The flux of the 0.1/im grain 
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Fig. 4. — Dust absorption cross sections for the three grain radii 
used in the exact radiative transfer calculations shown in figure 131 



TABLE 1 

The exact single-grain DDN model results 



Grain radius 




o{T.*)/(r(T sub ) 


i> 


q(2fim) 


gray 


49 


1.0 


4.0 


1.00 


2fim 


52 


1.3 


3.6 


0.99 


0.5^tm 


68 


3.6 


2.3 


0.45 


O.lfim 


150 


28 


1.3 


0.10 



model, which is almost angle independent as predicted 
by equation IA10I can reproduce only the weakest near- 
IR bumps. Therefore the DDN model can explain the 
near IR bump in Herbig Ae stars only when both of the 
following conditions are met: 

1) the rim dust opacity is gray in the near IR (grain 
radius > 0.5/im), 

and 

2) the disk is puffed to a height H r i m / ' R r i m ^ 0.15. 
Figure H shows that for these conditions to be satisfied, 
the DDN model requires the complete absence of small 
grains in the disk inner region. Therefore, for this model 
to work, the rim dust must undergo substantial growth 
that also fully depletes the population of small grain. At 
the same time, this process cannot be so extreme in the 
rest of the disk because the mid IR spectrum of Herbig 
Ae stars displays t he dust features of small grain emission 
l)van Boekel et al.ll2005j) . 

3. THE NEAR IR BUMP AND IMAGING EXPLAINED WITH 
A DUSTY HALO 

A dusty halo around the disk inner regions (~ 10 
AU) h as been invok ed to explain polarimetric measure- 
ments l|Yudinir2000() and correl ations between va riabili- 
ties in the optical and near IR (Eiro a et alJl2002j) . Such 
small regions are not yet accessible to direct imaging 
but have been resolved in near I R int erferometry by 
iMillan-Gabet. Schloerb. fc Traub I l|2001fl who also favor 
the halo geometry, although the interpretation of these 
visibility data is still model-dependent. Direct imaging 
is currently available only for larger scales, and these ob- 
servations have revealed larger halos, > 100 AU, around 
some objects (V03). The relation between the inner and 
outer halos, whether they are simply the inner and outer 
regions of the same circumstellar component, remains an 
open question. However, at the phenomcnological level 
this issue is not relevant because the two can be treated 



as separate circumstellar components if both are opti- 
cally thin. The inner halo is then radiatively decoupled 
from the cooler outer halo, simplifying the study of inner 
halos. 

Here we explore the contribution of the inner halo to 
the near IR emission. The halo precise geometry is not 
particularly important. It could be elongated, clumpy 
or inhomogeneous, but as long as it is optically thin 
it can be approximated with spherical geometry. The 
reason is that the temperature of optically thin dust is 
dominated by the stellar heating, resulting in spherically 
symmetric isotherms and circularly symmetric images 
at wavelengths where the dust thermal emission dom- 
inates over scattering (V03). Optically thin halos are 
also transparent to the disk emission and we can ignore 
the disk effect on the halo. The exact image shape ulti- 
mately depends on detailed dust density and grain prop- 
erties, telescope resolution and sensitivity, observational 
wavelength and the intrinsic ratio between the disk and 
halo surface brightness. Various observations of R M on 
vividly illustrate these effects (see lWeieelt et aLl l2002'). 

If the halo optical depth at visual wavelengths ry is 
larger than jH/R, where H/R is the disk flaring at the 
halo outer radius, then the halo dominates the SED com- 
ing from the dust within radius R around the star (see 
V03 for details). At near IR wavelengths, this condition 
is satisfied for the halo optical depths of interest here 
(Tv > 0.1). 

3.1. Theoretical Examples 

Our models consist of a star surrounded by a spherical 
halo with radial density profile r\ oc r~ p . The halo ex- 
tends from inner radius Rt n , set by the dust sublimation 
temperature T su b, to outer radius R ou t- 

The dust c hemist ry is x=0A olivine from 
iDorschner et a l. (1995), with grain size distribution 
n(a)oc a~ q between the minimum grain radius a m i n 
and maximum a max . We use q = 2, a m i n = 0.01/im 
and vary a max . The radial optical depth of the halo 
is specified at A = 0.55/^m as jy. The radiative 
transfer problem is solved with the code DUSTY 
l)Ivezic. Nenkova. fc ElitzuTI Il999t) . which takes advan- 
tage of the scaling properties of the radiative transfer 
problem for dust abs orption, emission and scattering 
i|Ivezic & Elitzurl ll997tL 

Figures [S] and [5] show some SED examples for dusty 
halos around 10,500K and 5,000K stars, representative 
of Herbig Ae and T Tauri stars, respectively. The stel- 
lar spectrum is taken from Kurucz models. In addition 
to the strength parameter F^/im/^i/xm, the flux ratio 
Fa^ito./ 'F2[im can be used to characterize the NIR bump 
shape. Both the strength and shape parameters are in- 
fluenced by changes in the dust sublimation temperature, 
maximum grain size, halo outer radius and optical depth. 
Comparison of the data with halo model results for the 
strength and shape parameters is shown in figure for 
the same objects as in figure Models for p = 1 and 
p = 2 halos around a 10,500K star are dispersed all over 
the diagram. Arrows show how the model results move in 
the diagram as the model parameters are varied, indicat- 
ing that various degeneracies are possible. The observed 
levels of bump strength and shape are readily reproduced 
with plausible values of the model parameters. 

We briefly summarize the effect of various halo param- 
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Fig. 5. — SED variation with the parameters of small spherical halos around stars with 10,500 K and 5,000 K (Kurucz stellar models). 
The radial density profile is r) oc y~ 1 . The dust chemistry is x = 0.4 olivine from Dorschncr ct al. (1995). The grain size distribution is 
n(a)oc a~ 2 between minimum grain radius 0.01/xm and maximum a max , as marked. The other varied parameters are the dust sublimation 
temperature T su t,, which sets the halo inner radius Ri„, the halo outer radius Rout and its optical depth at visual t. Note the variation of 
Ri„ among models, especially when a ma x is increased. 



Vinkovic et al. 



1 

0.1 
0.01 
10~ 3 
10~ 4 
10~ 5 

0.1 
0.01 



1 o -3 y 



10" 4 r 



10~ 5 b- 




10" 5 k- 



- 1 1 1 — I — I I I I 



1 1 1 1 — I — I I I I I 



- 1 1 1 — I — l l l I; 



1 1 I I I I I I | 1 1 I I I I I I | 1 1 I I I I I I : 



t=0.1. R„„,= 10R, = 999R. 



t=0.3, R out =10Ri 

_ _ r=0.6. R out =10R i 

T, ub =1100K 

Q moj< = 1/^ m 

10500K star 



-1030R, 
=1080R, 



H 1 I I I I 1 1 1 1 1 I I I I I 1 1 




T Iub =1100K, °.„= 1 i" m . R <, u t= 10R i„ = 1030R . 

T sub =1100K. a mo „=10/im, R oul =10 R in = 502R, 

_ T sub =1500K ' ° m „= 'l" m . R O ut =10R in= 438R . \ ■ 

T sub =1500K, a mo „=10/im, R oul =10R, n = 260R, 

1 1 I I I I I I | 1 1 I I I I I I | 1 1 I I I I I I : 




X (yt4m) 

Fig. 6. — The same as figure IS1 except that r) oc y~ 2 . 



100 



1 

0.1 

0.01 
10~ 3 

10- 4 

10~ 5 

0.1 
0.01 
10~ 3 

10- 4 

10~ 5 



- 1 1 1 — I — I I I I I 1 1 1 — I — I I I I I 

\ 



- 1 1 1 1 1 I I I; 



_ ° = 



t=0.3 
5000K star 



1 fim, 


R out= 


10.0R 


n =923R, 




R out= 


13.0R 


n=923R, 


5/im, 


R out = 


16.5R 


1 =923R, 


1 


R out= 


18.3R 


„=923R, 



"H 1 I I I I I I | 1 1 I I I I I I | 1 1 I I I I I I ; 



t'\\ 



t=0.1, R olll =10R in =888R, ' '\ 

T = 0.3, R„ 11( =10R,„ = 923R, •> 



t=0.6, R out =10R In = 973R, 

T sub =1100K 

□ maK = 1/zm 
5000K star 



A 

A 

'A. 



-H 1 I I I I I I | 1 1 I I I I I I | 1 1 I I I I I I ; 




X Gum) 



100 



Optically thin halos ari 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 



F A (4 M m)/F A (2 M m) 

Fig. 7. — Diagram of the strength (vertical axis) vs. shape (hor- 
izontal axis) of the near IR bump. Data are for the same sources 
as in figure Objects with enough data to determine the bump 
shape are marked with stars, otherwise only their bump strengths 
are marked with arrows on the right. Theoretical results for the 
halo models presented in figures 1ST and 151 are marked with crosses 
for halos around a 10,500K star and with T for halos around a 
5,000K star (the location of this naked star is marked with en- 
circled T). Model results for the 10,500K case are connected with 
arrowed lines that indicate their path in the diagram when only 
one parameter is varied, as marked. Thick lines correspond to p = 
2 halos, thin lines to p = 1. 



eters on the strength and shape of the near IR bump. 

Optical depth: A larger optical depth results in a 
stronger near IR bump. This reflects the dependence 
of flux on the total mass of emitting dust (equations A7 
and A12 in V03). The dust sublimation radius Ri n is 
only slightly affected, as expected in the optically thin 
limit where the diffuse radiation is negligible. 

Grain size: Larger grains shift the near IR bump to- 
ward longer wavelengths and make it appear more flat. 
With increased grain size the opacity becomes more sim- 
ilar to gray dust, resulting in a r~ ' 5 temperature profile 
since the geometrical dilution of stellar heating is the 
only cause of temperature variation. Smaller grain sizes 
create steeper radial temperature profiles. Therefore, for 
a given density profile smaller grains emit relatively more 
radiation at shorter wavelengths than larger grains. In 
practice, grain sizes come in mixtures and sublimate at 
different radial distances, greatly adding to the complex- 
ity of the problem. The SED models are therefore prone 
to various model degeneracies. 

Sublimation temperature: With a higher dust sublima- 
tion temperature, the near IR bump shifts to shorter 
wavelengths, reflecting the shift of the emission peak. 

Outer radius: The halo size can affect the near IR 
bump in two ways. On one hand, reducing the outer ra- 
dius while keeping the dust distribution fixed reduces also 
the total optical depth. The near IR bump then starts to 
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Fig. 8.— Models for the near IR bumps in AB Aur, HD 100546 
and HD 163296. Points and thin solid lines are the data (see 
Measurement uncertainties in the near IR are comparable to the 
symbol sizes. Thick solid lines are the models, comprised of the 
stellar component (dotted lines) and a halo contribution (dash- 
dotted lines); see text and table|z]for model details. In the case of 
AB Aur, the thick dashed line is the second model listed in table |21 
Only data at wavelengths shorter then 6^tm were employed in the 
fits since in this study we consider only the inner halo responsible 
for the near IR bump. Longer wavelengths (the gray area) are 
displayed to show the halo contribution to the mid IR features; in 
HD 163296, this contribution suffices to explain the observed 10/jm 
feature. 



decrease when the dust removal reaches the NIR emis- 
sion regions at radial distance ~ 10i?i„ (temperatures > 
500 K). On the other hand, reducing the outer radius 
at a fixed halo optical depth is equivalent to redistribut- 
ing the dust within the halo. The bump then becomes 
stronger as the outer radius is reduced because more dust 
is shifted toward smaller radii and higher temperatures. 

Stellar temperature: As its temperature decreases, the 
emission from the star starts to blend with that from the 
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Fig. 9. — Sketch of the circumstcllar geomet r y indi c ated by imag- 
ing ob serva tions o f HD 100546 IGradv et aT] 120011) . IGradv et all 
(2005); see H3.2.2I . Region A is an optically thin dust layer, dom- 
inating the near IR flux. The disk, marked with B, is cooler and 
does not affect the near IR wavelengths. Both regions A and B are 
embedded in an optically thin dusty envelope, marked with check- 
ered pattern, whose optical depth is so small that its contribution, 
too, to the near IR flux is negligible. 



halo, and the near IR bump disappears. Only a careful 
analysis can then separate the stellar from the diffuse 
flux in the near IR and reveal the bump. For comparison 
with T Tau stars, figure shows also 5,000 K models 
(marked with T). In spite of the large variations in halo 
parameters, these models display only a limited range of 
bump strengths and shapes close to the naked star values. 
This explains why the near IR bump was not originally 
recognized in T Tau stars while easily detected in Herbig 
Ae/Be stars. 

3.2. Observational Examples 

We show three examples which illustrate different cir- 
cumstellar dust configurations: AB Aur, HD 100546 and 
HD 163296. HST imaging suggests that AB Aur and HD 
100546 have large halos at radii > 1 " , while HD 163296 
shows only a disk ijGradv et alJ20 03) . Irrespective of the 
existence of a large halo, all three objects show a near IR 
bump, with the strongest bump in AB Aur. Since the 
focus of this study is the near IR bump, the large scale 
halos are irrelevant here and we only consider a small 
halo within ~10 AU around the star. 

Our fits to the data are shown in figure |H1 with the 
model parameters listed in table [21 The halo outer ra- 
dius is 10 times the dust sublimation radius in all models. 
Since our focus is the near IR bump, our model consists 
only of the star and the inner halo, and only data at 
wavelengths shorter then 6/xm were employed in the fit- 
ting. The derived model parameters are not unique since 
various degeneracies exist in model results for the near 
IR flux (sec < j3.1l) . For exa mple, the "hot component" 
in the lBouwman et"al] l)2000|) models can be interpreted 
as a small-scale halo with dust properties different from 
those in our study. 

3.2.1. AB Aur 

The emission from AB Aur has been resolved at var- 
ious wavelengths and interpreted as a disk with vastly 
different estimates for the inclination angle, as follows: 



visual: i < 45° Oradv et al. (19991 

near IR: i < 30° Eisner et al. (20031 . (2004) 

i = 30 ± 5° Fukagawa et al. (20041 

mid IR: i = 55 ± 10° Liu et al. (20051 

millimeter: i = 17°tt Semenov et al. (2004) 
i = 21°.5tu-3 Onrder et al. (20051 
% = 33 ± 10° Pietu et al. (20051 
i ~ 76° Mannings fc Sargent (19971 

Such a disparity is expected in halo-embedded 
disks (see figure 7 in V0 3) because, as noted by 
iMiroshnichenko et alJ l)1999D . the halo dominates the im- 
ages at wavelengths extending to ~ 100 /im or so, and the 
disk emerges only at longer wavelengths. Interpretation 
of molecular line images, too, must be done carefully to 
avoid confusion between the halo and disk contributions. 

A general conclusion about the AB Aur inner halo is 
that it must have a radial density profile between 1 /r and 
1/r 2 ; this differs from the outer halo, which has a 1/r 
density profile as deduced from the 1/r 2 radial b right- 
ness profile of the HST image llOradv et al.l l"l999. sec 
also equation A10 in V03). Conclusions regarding the 
properties of the dust grains in the inner halo are less 
firm. Near IR interferometry suggests the presence of 
dust close to the star, implying large grains that can sur- 
vive at small distances. An example of a big grains model 
for AB Aur is shown in figure OH with thick dashed line 
(see also table |2J. The grain size and chemistry might be 
subject to radial variati ons, as is indicated by compar- 
ison bet ween the HST llOradv et alJ 119991) and Subaru 
images l|Fukagawa et alJl2004D .' further complicating the 
modeling. 

3.2.2. HD 100546 

The HST image of this source l|Gradv et all l200l|) 
shows a very tenuous large scale nebulosity, whose low 
surface brightness implies an optical depth of only ry ~ 
0.015. This component of the dust distribution does 
not contribute significantly to the IR emission and can 
be ignored in the current analysis. The HST image, 
which is produced purely by scattered light, reveals also 
a p rominent disk with in clination angl e 49 ± 4°. Near 
IR ijAugereau et al.ll200lD and mid IR ijLiu et alJl2003T) 
imaging give similar results for the disk even though the 
latter is produced purely by dust emission and the for- 
mer contains a mixture of both scattering and emission. 
The HST brightness contours are symmetric, with the 
brightness declining as 1/r 3 . These are the signatures 
of scattering from the CG layer of a flat disk (see V03). 
However, for the CG model to be applicable, every point 
on the scattering surface, which extends to a distance of 
~10" from the star, must have a direct line of sight to 
the stellar surface. This is impossible in the case of a flat 
disk, since it would have to maintain a thickness smaller 
than the stellar radius for hundreds of AU. Therefore, 
the only self-consistent explanation of the HST imaging 
is with an optically thin halo whose dimensions are un- 
related to the stellar size. The HST image implies that 
the halo has a flattened geometrical shape, and its 1/r 3 
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TABLE 2 

Halo model parameters for three case studies of Herbig Ae stars 



Object 


T* 


Av 


T sub 


TV 


V 






Fbol 


grain radius 


carbon 


olivine 


AB Aur 


9750 


0.35 


1500 


0.35 


r~ 


2 


0.78 


0.9 


q = 3.5, a max = 0.25 


40% 


60% 


AB Aur 


9750 


0.35 


1800 


0.35 


r~ 


2 


0.25 


0.9 


q = 3.5, a max = 5.0 


40% 


60% 


HD 100546 


10500 


0.3 


1500 


0.35 a 


r~ 


1 


0.45 


0.8 


q = 2, amax = 0.50 


10% 


90% 


HD 163296 


9500 


0.3 


1500 


0.25 


r~ 


1 


0.61 


0.9 


q — 3.5, CLrnax = 0.25 


30% 


70% 



Note. — Description of columns: T* — stellar temperature in K. T su b — dust sublimation temperature in 
K. Av — reddening toward the star by dust other than the inner halo. Tv — visual optical depth of the halo, n 
— radial dust density profile of the halo. Ri n — halo inner radius in AU, determined from dust sublimation; the 
halo outer radius is lQRi n . Fbol — total bolomctric flux in 10 ~ 1D W/m 2 . grain radius — dust size distribution 
a~ q between minimum q m j n — 0.005 ftm and maximum a max listed in )im. Amorphous carbon properties from 
IHannerl I1988D . olivine from lDorschner et afl 119951) (with x = 0.4). 

a HD 100546 is modeled with a flattened halo, which docs not contribute to the circumstcllar reddening because 
the dust is out of the line of sight toward the star (see figurclol . 



brightness profile implies that it has a 1 /r 2 radial density 
profile (V03). This flattened halo is outlined as region 
A in figure El Since the halo dominates the imaging, the 
geometry of the optically thick disk structure, outlined 
as region B in the figure, remains unknown. 

The HST imaging does not constrain the inner-halo 
geometry at radii < 10 AU. The surface density must 
be reduced in that region because the near IR bump in 
HD 100546 is significantly smaller then in AB Aur even 
though otherwise the two stars are rather similar. In- 
deed, the fit to the near IR bump yields al/r radial den- 
sity profile (figure[H]and tableEJ) , shallower than in the re- 
gion resolved by HST. The fit was further improved by an 
increased contribution from large grains (a max — 0.5^m 
and q = 2) and a redu ced fraction o f carbo n dust in the 
mix. Observations by iGradv et alJ l)2005[) show that a 
constant density profile, creating 1/r brightness profile, 
might be more appropriate in the region between 20 and 
50 AU. We find that a constant density model could also 
fit the near IR spectrum if the sublimation temperature 
were increased to 1700 K. All these results point toward 
large structural differences between the inner and outer 
regions of HD 100546. 

3.2.3. HD 163296 

The model properties of the inner halo in this source 
are very similar to AB Aur, except that a shallower den- 
sity profile of p — 1 is preferred (see table [5J ■ A similar 
general conclusion is that the halo radial density profile 
is between 1/r and 1/r 2 , with uncertainties in the dust 
properties. Significantly, in this source the inner halo 
also fits the 10/im feature all by itself (figure [BJ. No 
other optically thin components are required for explain- 
ing the mid IR dust features, indeed none are observed. 
The HST image, which is incapable of resolving the inner 
halo, shows no evidence of a large scale str ucture other 
then the disk, with an inclination of 60 ±5° l|Gradv et alJ 
2000) in agreement with ~ 58° derived from millimeter 
observation l)Mannings fc Sargenllll997lh Another note- 
worthy feature of the HST image is a bipolar jet. The 
process responsible for jet formation could perhaps also 
lift up dust above the disk and create the small scale 
halo responsible for the near IR bump. Such possible 
correlation between jets and the near IR bump can be 
studied further when more high resolution data from a 
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Fig. 10. — Size-luminosity diagram for the inner-halo model. 
The sublimation radius Ri n of each displayed object is determined 
from a simultaneous fit to its near-IR bump and visibility data. 
Filled circles mark objects with dust sublimation temperature 1500 
K, empty circles with 1800 K. Lower limits on Ri n are shown with 
solid line for T auB = 1500 K and dashed line for T auo = 1800 
K; they correspond to halos of gray dust and zero optical depth. 
The dotted line connects a group of objects that follow closely the 
1 /2 

Ri n oc LJ relation, indicating that their inner halos have similar 
optical depth and dust properties. The original version of this 
diagram, constructed from rin g-model visibility analysis, prod uced 
a pu zzling amount of scatter (Monnicr & Millan-Gabct 2002, see 
H3.3I . The inner-halo model resolves the puzzle. 



larger sample of objects become available. 

3.3. Size-luminosity Correlation 

The milli-arcsccond resolution reached in near IR inter- 
ferometry enables studies of the immedia te env ironment 
of young stars, down to 0.1 AU l)Malbetll2003|) . Unfor- 
tunately, current visibility data are not yet capable of 
reproducing the full 2D image of an object, instead re- 
quiring a model of the geometry for their analysis. One 
simple and often used model of the circumstellar geome- 
try is a flat dust ring of uniform surface brightness. This 
ad hoc model did not arise from some specific radiative 
transfer modeling but rather chosen as a simple approach 
to the visibility fitting procedure. 

Fitting the visibi lity data of a number of obje cts with 
this ring model, iMonnier fc Millan-Gabetl l)2002|) discov- 
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ered that the size of the ring inner radius increased with 
the stellar luminosity L*. This is the expected result 
when dust sublimation controls the size of the dust-free 
region aro und the star. Since rad iative transfer is scale 
invariant l(Ivezic & Elitzurl 119971) , inner radii of rings 

1/2 

would be expected to scale as LJ if their dust proper- 
ties w ere the same. However, iMonnier fc Millan-Gabetl 
( 2002) do not find such a trend. Instead, at a fixed lu- 
minosity the derived radii vary by almost a factor of ten, 
which they refer to as scatter in the size-luminosity di- 
agram. This scatter indicates either that the disk inner 
regions have vastly different properties, with the sub- 
limation temperature varying from ^1000 K to ^2000 
K, or that the ring model is n ot a proper description of 
the ac tual dust distribution. IMonnier fc Millan-Gabetl 
(2002) also noted that some highly luminous objects 
(L* > 10 3 Lq) had smaller than expected inner ring radii, 
thus requiring even higher dust sublimation tempera- 
tures. New interferometric data bv IMonnier itjjl] (2005) 
slightly reduce the scatter in the ring-radius-luminosity 
relation, but the remaining scatter still implies a large 
range of sublimation temperatures, and very luminous 
objects still display abnormally small radii. 

Instead of the ring model we have analyzed the inter- 
ferometry results with the inner-halo model, performing 
simultaneous fits of both the near-IR bump and visibil- 
ity data. Preliminary results are shown in figure IT771 It 
is highly significant that there are no objects in the for- 
bidden region below the indicated lower limits. The cor- 
relation of overall bolometric luminosity with inner ra- 
dius is much tighter than in the ring model, the small 
remaining scatter arises from variation in halo optical 
depth and grain size. In contrast with the ring model, 
the sublimation temperature rarely differs from 1500K 
(it is 1800K in a couple of objects). The high luminosity 
object MWC297 which was esp ecially troubling in the 
IMonnier fc Millan-Gabetl 1)20021) analysis is now consis- 
tent with 1500K sublimation temperature. It is striking 
how some of t he objects that were highly sca ttered in 
the diagram bv lMonnier & Mill an-Gabetl (2002) now set- 

1 /2 

tie on the same Lj size-luminosity relation (dotted line 
in figure I10|l , indicating similarities in the halo proper- 
ties of all these stars, in turn pointing toward a common 
physical mechanism of halo formation. 

4. CONCLUSIONS 

An examination of the puffcd-up disk rim model 
(DDN) shows that it has rather limited capabilities in 
explaining the near IR bump of Herbig Ae/Be stars. The 
observed level of near IR excess implies a certain emit- 
ting volume of optically thin puffed-up disk rim surface 
for given dust properties. The volume derived from the 
DDN model falls short of this observational limit, unless 
the disk is made of perfectly gray dust. The puffed-up 
rim produces enough near IR flux only when the inner 
disk consists purely of dust grains larger than ~ 0.5 /im 
and the disk puffin g reaches values of H r j m I R r im ^ 0.15. 
Models bv lDullemond fc Domimkl l|2004l) show that even 
traces of small grains inhibit the disk puffing, eliminat- 
ing the DDN model as a viable explanation of the near 
IR bump. Since the 10/xm emission feature indicates the 
presence of small grains in the circumstellar dust, addi- 
tional mechanisms must be invoked to remove all small 



grains from the inner disk and keep the DDN model vi- 
able. 

From fits to the SED of a number of Haebes, 
iDominik et alJ l)2003j) conclude that the infrared excess 
in these stars is produced by disks alone without the 
need for additional circumstellar components. This con- 
clusion is invalidated by the mathematical proof that a 
fit to the SED cannot distinguish between the surface of 
a flared disk and an optically thin halo (V03). Fits to 
the SED alone are not a conclusive proof of a particular 
dust geometry. 

We find that that the optically thin dusty halos around 
the disk inner regions whose existence has been inferred 
in various observations readily explain the strength and 
shape of the near IR bump. The halo is not limited by the 
disk properties. Hence, it can extend above the disk sur- 
face and accommodate the emitting optically thin dust 
volume required by the near IR flux observations. The 
required halo is rather small, less than several AU in size, 
and its optical depth in visual is less than ~ 0.4. Despite 
its small optical depth, the halo dominates the near IR 
spectrum and hides the disk near IR signature. However, 
detailed properties of the halo, such as its exact shape, 
grain properties or dust density profile, are not uniquely 
constrained by the SED since different combinations of 
the parameters can produce the same flux. These de- 
generacies can be broken only with imaging capable of 
resolving the disk inner regions. 

Inner halos not only explain the near IR bump but 
also successfully resolve the puzzle presented by the re- 
lations b etween luminosities and near-I R interferomet- 
ric sizes ijMonnier fc Millan-Gabet1 l2002). In addition to 
their near IR emission, the halos contribute also to the 
mid IR flux. HD 163296 is an extreme example where the 
halo in itself fully explains the mid IR dust features with- 
out the need for additional extended components (figure 
IHJ . The absence of such components in the HST image of 
this source is another success of the inner-halo model. In 
general, though, the inner halo emission is not expected 
to dominate the mid IR but still make a significant con- 
tribution that must be included in fits to the overall SED 
for reliable mod elling of the rest o f the ci rcumstellar ma- 
terial. Recently Ivan Boekel et, all il200M ■ suggested that 
differences in the strength and shape of the mid IR sil- 
icate feature in Haebes are evidence for dust settling in 
the disk. However, these differences could instead reflect 
halo evolution, with the most active stars showing the 
strongest mid IR signature of the inner disk halo. High 
resolution imaging is necessary for definite conclusions 
about evolution either of the dust or the circumstellar 
disk. Such imaging will soon become available at the 
VLTI, which offers milli-arcsecond resolution at near IR. 
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APPENDIX 

A. APPROXIMATE SOLUTION FOR THE RIM RADIUS 

The observed emission from a puffed up inner disk depends on the rim radius, height and surface brightnes s (see |)2.2I). 
Derivi ng the rim radius requires a proper treatment of the rim temperature structure. As already noted bv lDullemondl 
( 2002) , in gray dust the diffuse radiation creates a temperature inversion — the dust temperature is maximum in the 
rim interior (that is, at R > R r im), riot on the rim surface. Here we derive an approximate solution for the dust 
temperature T on the rim surface and the temperature T\ inside the rim at depth Ty ~ 1 from the surface. The 
solution demonstrates the inversion effect for gray dust and shows that it does not exist in the non-gray case. 

A.l. Dust temperature Tq on the rim surface 

Consider a dust grain on the rim surface. It is heated by the stellar flux and the diffuse flux F out coming out 
from the rim interior. For large optical depths these two fluxes are balanced: F* = F out . Stellar flux absorbed by the 
grain is ayi 7 "*, where ay is the dust cross section in visual. Absorbed diffuse flux is ~ 2ai R F out , where ai R is the cross 
section in the near IR and the factor 2 accounts for absorption from 2ir steradian. The grain emits into Air steradian, 
so that the energy balance is 

ayF* + 2a IR F out = Aai R a S BTo (Al) 
where asB is the Stefan-Boltzmann constant. Using F» = F out we get 

a SB T* = E±( 2 + ^L) (A2 ) 

A. 2. Dust temperature T± at Ty ~ 1 from the surface 

Now consider a dust grain at distance Ty ~ 1 from the surface into the rim. This grain is heated by the attenuated 
stellar flux F* exp(— Ty) = F„ exp(— 1), by the diffuse flux from the surface dust between Ty = and Ty ~ 1, and by 
the diffuse flux from the rim interior. The absorbed stellar flux is ayF* exp(— 1). Diffuse contribution from the surface 
dust layer is ~ 2a i R a sbTqTir-, where tir ~ am/ay is the infrared optical depth of this surface layer. Diffuse heating 
from the rim interior, described by temperature Ti, is ~ 2aj R as B Ti- This is a good approximation for gray dust and 
an overestimate for non-gray dust, where that temperature decreases rapidly with optical depth. The energy balance 
is 

a v F* exp(-l) + 2a IR a SB TQT IR + 2a IR a SB T? = Aai R a SB T{ (A3) 
Using To from equation IA2I we get the interior temperature 

* aB Tt=*±Um + \ + 2*LeA (A4) 
4 V °v ctir J 

Note that the ratio Tq/T\ depends only on ai R /ay and is independent of F*. 

A. 3. Gray and non-gray regimes 

We consider two distinct opacity regimes: gray when ai R /ay ~ 1 and non-gray when ai R /ay <C 1. The ratio of the 
rim surface temperature To and the interior temperature T\ in these two regimes is 

Tq/Tx = 0.95 when a IR /a v = 1 (gray dust) (A5) 
T /Ti ~ 1.08 when a IR /a v < 1 (non-gray dust) (A6) 

The gray opacity creates a temperature inversion with the temperature in the rim interior higher than on its surface. 
This inversion does not appear in non-gray dust where the temperature decreases monotonically with distance from 
the rim surface. If the maximum dust temperature is 1,500K (sublimation temperature) then gray dust has Ti = 1,500 
K and T -1,400 K, while non-gray dust has T =l,500 K and T x <1,400 K. 

These approximate expressions are in reasonable agreement with the results of exact 2D radiative transfer calculations 
(see H2.3|l . which yield 

T =1387 K and Ti = 1490 K for gray dust 

T o =1500 K and Ti = 1236 K for O.l^m grains 

The transition between these two regimes occurs at grain radius 0.5/Ltm which yields To=1474 K, Ti = 1492 K and the 
maximum temperature of 1,500 K at Ty =0.34. 
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A. 4. Disk rim radius and near-IR bump strength 

Since the rim optical depth is large, we can assume that both temperatures To and T\ are located at essentially the 
same distance from the star. If we set T\ to the dust sublimation temperature T su b then based on equation IA4I the 
rim radius is 



R rim = -Rj^) */2— + 1 4- 2-^e- Vl + H„ m /R rtm (A7) 

& X^subJ V °V &IR 



where we used T» = L* ^1 + H r i rn / R T i m / 4:7rRf im (the factor y/T+ HrimjRrim is correction for rim self- irradiation, 
already introduced by DDN). Comparison with equation 0] for the rim radius gives 



^2 + ^fl + 2^V (A8) 
e cry V OV / 



The two extreme opacity regimes yield 



{1 when am/ cry <C 1 (non-gray dust) 
(A9) 
4 when (Jir/<Jv —* 1 (gray dust) 

This result can also be derived by setting To ~ Tsub- Combining this result with equations [31 and 01 and dividing the 
overall observed flux at 2/xm by the stellar flux at 1/xm yields the near-IR bump strength 

TP R IT \ / T 1 \ 4 R It \ ff / ff \ f 1 non-gray dust 

r 2[±m JD 2l^m\- L * ) , / ^ * \ J ^2p 1 Tn \± rim) ^rira I ^ . ^-rira \ 1 

^lpm -^l^m(T*) \Ty^ m y -B^^ im (T t: ) 7rR r j yTn y Rrim ) 



(A10) 

4 sin i gray dust 



where we used the approximation c/>, ~ <7ir/<7\/ ~ cr (T r i m )/<r(T (! ). This solution shows that non-gray dust gives 
angle-independent bump strength in addition to reducing its magnitude from the gray dust result. 
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